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We theoretically investigate optical injection in semiconductor ring lasers and disclose several 
dynamical regimes. Through numerical simulations and bifurcation continuation, two separate 
parameter regions in which two different injection-locked solutions coexist are revealed, in addition 
to a region in which a frequency-locked limit cycle coexists with an injection-locked solution. Finally, 
an anti-phase chaotic regime without the involvement of any carrier dynamics is revealed. Parallels 
are drawn with the onset of chaos in the periodically forced Duffing oscillator. 
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I. INTRODUCTION 

Optically injected laser systems generally consist of 
two laser sources, a "master" laser whose output light is 
coupled into the cavity of a second "slave" laser. A sim- 
ple model for this type of system is a nonlinear oscillator 
(slave) which is periodically driven (master). Although 
these systems are relatively simple, they exhibit a wealth 
of dynamical behavior which has been widely studied for 
different types of lasers [TrUO] . 

A class of semiconductor lasers for which the nonlinear 
dynamics induced by optical injection have not yet been 
investigated in depth are the semiconductor ring lasers 
(SRLs). A SRL is a semiconductor laser in which the 
light is confined in a circular waveguide structure. As 
a result, SRLs generate light in two opposite directions 
referred to as the clockwise (CW) and the counterclock- 
wise (CCW) mode (see Figure [T]) . SRLs have received 
increasing attention in recent years because they 

are suitable candidates as key components in photonic 
integrated circuits [T2]. The bistable character of their 
directional mode operation allows them to be used in sys- 
tems for all-optical switching and as all-optical memories 
|12[ I13j . Optical injection can be particularly important 
in SRLs as a control mechanism for the dynamics in op- 
tical switching applications [21 [TS] , or when a holding 
beam is used to enforce unidirectional operation in the 
injected direction |16| . However, optical injection can 
also give rise to very intricate dynamics, which may ob- 
struct the desired dynamical behavior. The particular 
Z2-symmetry of the SRL and its resulting phase space 
structure has already led to results such as alternative 
switching mechanisms [T51 [T7], multistable regimes [TS] 
and non-Arrhenius mode-hopping |19j . The different dy- 
namical regimes resulting from forcing the SRL through 
optical injection will be disclosed in this paper. 

Given the interest in SRLs, we will investigate the be- 
havior of this two-mode device when subjected to uni- 
directional optical injection. Optical injection in other 



two-mode semiconductor lasers such as VCSELs [H] and 
two-color lasers [TU] has also recently been investigated. 
The bimodal character of these devices gives rise to more 
intricate dynamical behavior, such as bistable operating 
regimes governed by optical injection. 

In this paper we will perform an extensive bifurcation 
analysis using the software continuation package AUTO 
[20 , complemented by numerically solving the rate equa- 
tion model. This approach allows us to compute the bi- 
furcation diagrams for optically injected SRLs and to re- 
veal the stability of the invariant structures in the phase 
space of the mathematical model which will be intro- 
duced in the next section. Keeping in mind the all-optical 
memory and switching applications, we focus on unidi- 
rectional optical injection. This means we optically inject 
only one of the two counterpropagating modes, chosen to 
be the CW mode. This will also facilitate comparison of 
our results with future experiments. We will further as- 
sume that the SRL is biased in the bistable unidirectional 
regime. 

This paper is organized as follows. The rate equation 
model of the optically injected SRL is described in Sec- 
tion [TTJ where we also point out the differences between 
the SRL model and that of other two-mode lasers. The 
analysis starts in Section |III| in which we reveal some of 
the characteristic behavior of the optically injected SRL 
obtained by numerically solving the rate equation model. 
In Section |IV[ we complement this analysis by present- 
ing the bifurcation curves of stationary points, and point 
out the differences with other optically injected semicon- 
ductor lasers. This allows us to give a bird's eye view 
of the dynamics exhibited by the optically injected SRL 
for different values of the detuning and injection powers. 
For a certain parameter range, the optically injected SRL 
exhibits a novel anti-phase chaotic regime, described in 
Section [Vj We finally draw conclusions of our analysis 



and point to future work in Section VI 



II. FORMULATION OF THE MODEL 
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We consider a typical master-slave setup in which we 
neglect coupling from the slave to the master laser (see 
Figure 111. In this setup, the SRL is assumed to operate 
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FIG. 1. (Color online) Schematic setup of the optically in- 
jected SRL (m is assumed to be the amplitude coupling coef- 
ficient of the straight output waveguide into the ring cavity) . 



in a single transverse and single longitudinal mode and 
can sustain two counterpropagating directional modes. 
Following Refs. [22 [22] with a straightforward modifi- 
cation to account for the optical injection as in [15], we 
can write the following rate equations for an optically in- 
jected SRL; neglecting spatial variations within the laser 
and adiabatically eliminating the medium's polarization 
dynamics: 
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Here t is time, E\ and E 2 are the slowly varying com- 
plex envelopes of the counterpropagating waves, TV is the 
carrier population inversion, [i is the renormalizcd injec- 
tion current (/i = at transparency and \i = 1 at lasing 
threshold), k is the field decay rate, 7 is the carrier decay 
rate, a is the linewidth enhancement factor and Ti n is the 
cavity round-trip time. The two control parameters are 
the injected field amplitude Ei > and its detuning A 
from the longitudinal mode frequency of the SRL. A de- 
tuning A > corresponds to a higher master than slave 
frequency. The two counterpropagating waves saturate 
both their own and each others gain through spectral 
hole burning and carrier heating effects. These self- and 
cross-saturation phenomena occur on faster time scales 
than the photon lifetime of the SRL [23 , allowing them 
to be added phenomenologically, modeled by s and c. 
Note that the cross-saturation is stronger than the self- 
saturation (c ~ 2s) |24j . In addition to this nonlinear 
coupling there also exists a linear coupling between the 
counterpropagating waves, referred to as backscattering. 
It is caused by reflections inside the cavity at the inter- 
face with the coupling waveguide and at the cleaved end 
facets of the output waveguide. They result in a linear 



coupling between the two fields, modeled by an ampli- 
tude k and a phase shift (pk- Finally, note that the ref- 
erence frame of the equations is chosen to corotate with 
the phase of the master laser so that fixed points of this 
system correspond to injection-locked states. 

In a typical experimental setup, the photon lifetime 
T p = k~ and the carrier lifetime r c = 7" 1 are respec- 
tively of the orders 10 ps and 5 ns, yielding two diffe- 
rent time scales in the system. The other parameters are 
fixed to realistic values a = 3.5, s = 0.005, c = 0.01, 
k = 0.4412 ns" 1 , cf) k = 1.4966 and r m = 0.6 ps [21]. The 
value of the bias current /i is chosen such that the SRL 
operates in the bistable unidirectional regime, but still 
relatively close to the alternate oscillation regime [17J. 
The detuning is varied up to 7 ns" 1 (angular frequency), 
while the values used for the injection amplitude Ei span 
several orders of magnitude, ranging from O(10 ) up to 

o(io- 2 ). 

For future reference, we will briefly highlight the char- 
acteristic timescales of the solitary SRL. The principal 
time scale in any semiconductor laser is the relaxation 
oscillation timescale. For our set of equations ([lj, the 
relaxation oscillation angular frequency can be approxi- 
mated by 



V^Ou - 1)jk w 5.307 ns" 



(2) 



In the case of SRLs there is also a second important 
time scale regarding intensity oscillations, the alternate 
oscillation frequency. It characterizes a particular op- 
erating regime of the solitary SRL in which two cou- 
pling mechanisms — the cross-gain saturation and the 
backscattering — compete with each other, inducing in- 
tensity oscillations at an angular frequency given by |21j 



loao = 2fc\/— cos(20 fc ) — asin(20 fc ) ps 0.606 



ns 



(3) 



In the bistable unidirectional regime, the SRL has four 
different steady state solutions in the absence of optical 
injection. Two of them are stable quasi-unidirectional so- 
lutions with the power concentrated in either the CW or 
the CCW mode. The other two are unstable bidirectional 
solutions with equal power in both modes, with the coun- 
terpropagating fields respectively in-phase (IP) and out- 
of-phase (OP). They are each characterized by a partic- 
ular optical frequency u x , with X={CW,CCW,IP,OP}. 
This is due to the different carrier densities associated to 
each of these solutions. The carrier densities are influ- 
enced by the optical intensity, in which the backscatter- 
ing also plays a role by altering the effective gain of both 
modes. This frequency corresponds to a certain detuning 
with respect to the cavity resonance frequency u which 
can be calculated numerically and is given by 



A cw = A ccw » 0.225 ns" 1 
A p = -Aip « 0.326 ns" 1 . 



(4) 
(5) 



with Ax = ujx — loq. Note that Acw must be equal to 
Accw due to the symmetry properties of our system (E\ 
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FIG. 3. (Color online) Representative time trace of the re- 
spective intensities of the CW (dark gray, red online) and the 
CCW mode (light gray, green online), and the total intensity 
(black) inside the cavity of the SRL when operating in the 
chaotic regime. Transient behavior has been eliminated, /i = 
1.704, fc = 1.4966, k = 0.4412 ns -1 , a = 3.5, c = 0.01, s = 
0.005, k = 100 mr\ 7 = 0.2 ns -1 , Ei = 6 10~~ 5 , A = -1 ns" 1 . 



FIG. 2. (Color online) Orbit diagram depicting the local ex- 
trema of the intensity |-Ei| 2 of the CW mode versus the injec- 
tion amplitude Ei. Black (gray) dots indicate local maxima 
(minima), /i = 1.704, 4> k = 1.4966, k = 0.4412 ns" 1 , a = 3.5, 
c = 0.01, s = 0.005, n = 100 ns" 1 , 7 = 0.2 ns" 1 , A = -1 ns" 1 . 



and E2 may be exchanged). The fact that Aop = — Aip 
follows immediately from the equations ([I]). 

Contrary to two-color lasers where the mode spac- 
ing is highly nondegenerate |10j . the mode spacing in 
SRLs is degenerate since the counterpropagating modes 
have identical frequencies. For this reason, the coun- 
terpropagating modes in SRLs have a significant phase 
coupling. On the other hand, the phases of the two 
modes in two-color lasers are decoupled and only the in- 
tensity of the uninjected mode influences the dynamics 
[TU] . Moreover, contrary to SRLs, saturation effects in 
two-color lasers are such that self-saturation is stronger 
than cross-saturation [10J . affecting the relative stabil- 
ity of the modes. Both the SRL model and the spin-flip 
model for VCSELs [251 account for a phase coupling 
between the modes. However, in the spin-flip model, the 
optical injection studied in Ref. [9] corresponds to injec- 
tion in both modes when the VCSEL is lasing in both 
circular modes simultaneously. Finally, the SRL model 
does not need an extra dimension for a second carrier 
population as the VCSEL model does. 



III. NUMERICAL SIMULATIONS 

We start our analysis by numerically solving (JlJ for 
different values of the detuning A and the injected field 
amplitude E i} using a fourth order Runge-Kutta algo- 
rithm with a fixed time step of 1 ps. We have constructed 
orbit diagrams for different values of the detuning using 
Ei as a parameter. These orbit diagrams plot the local 
extrema of the system's attractor as a function of the in- 
jection amplitude (or only a part of the attractor if its 



basin of attraction is not the whole phase space) . 

The orbit diagram for A = — 1 ns" 1 , only depicting the 
local extrema of the intensity of the CW mode, is shown 
in Figure [2j For low injection amplitudes the output 
power starts to oscillate around the former steady state, 
due to the beating between the injected signal and the 
optical fields inside the SRL. When increasing Ei, the 
amplitude of the oscillation grows and eventually the SRL 
switches mode (at point A) due to the optical injection 
in the CW mode (the solitary SRL was assumed to reside 
in the CCW mode prior to the optical injection). 

At point B the attractor changes from a limit cycle to 
a more complicated structure, yielding a sudden burst 
of local extrema. The time trace corresponding to this 
particular region is shown in Figure [3j In Section [V] we 
will later show that it corresponds to a chaotic regime. 
However, this regime has a different origin than those 
observed in other optically injected laser systems. From 
point B up to point C, the orbit diagram reveals chaos 
interspersed with periodic windows. 

For higher Ei (C— >T>), the system relaxes to a stable 
limit cycle which eventually dies out in a Hopf bifurcation 
at point D. At that point, the SRL locks to the injected 
signal. In the next section, we will reveal that this partic- 
ular injection-locked solution is only one out of three dif- 
ferent injection- locked solutions (the bidirectional one). 
The other two (unidirectional) injection-locked solutions 
have the same route to locking as the regular injection 
problem, which happens through a saddle-node bifurca- 
tion for low injection powers [21 HI [7]. 

Figure [2] shows a Hopf route to locking of the bidi- 
rectional injection-locked solution. There is also a Hopf 
route to locking in other optically injected semiconduc- 
tor lasers [3], but it only occurs for much higher injection 
powers. In that case, considering that for a fixed value of 
the detuning one would continuously raise the injection 
power from zero, the slave laser would first injection- lock 
through a saddle-node bifurcation, after which it unlocks 
because of the undamping of relaxation oscillations and 
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FIG. 4. (Color online) Bifurcation diagram in the (A, Em- 
plane, generated using the rate equations |TJ. Saddle- node 
(Hopf) bifurcations are depicted in black (light gray, red on- 
line). Supercritical (subcritical) bifurcations are depicted in 
full (dashed) lines. 



finally locks again through a Hopf bifurcation. In our 
scenario, the very first locking event happens through a 
Hopf bifurcation. 



IV. BIFURCATION ANALYSIS 

The dynamical behavior of the solutions of ([lj will 
generally vary for different values of the injected field 
amplitude Ei and the detuning A. Qualitative changes 
in the dynamical behavior of the system, so-called bi- 
furcations, can be numerically detected and continued in 
this two-dimensional parameter space using for example 
the bifurcation continuation package AUTO [20] • Figure 
[4] shows different regions in the (A,Sj)-plane bounded 
by bifurcation lines, with each region corresponding to 
different dynamical behavior. 

Two of these bifurcation lines, SNi and Hi, are famil- 
iar. They arise in many other optically injected lasers for 
small injection amplitude and detuning [5J [3J UJ (STJ 12H]- 
The SNi-line represents a saddle-node bifurcation on a 
cycle (infinite-period bifurcation), while the Hi-line rep- 
resents a Hopf bifurcation. The region confined between 
these two lines is the stable locking region, where the 
SRL is phase-locked to the injected signal, yielding an 
injection-locked solution s cw . The transformation of this 
solution near the boundaries of the stable locking re- 
gion is identical to other optically injected laser systems 
[2H]- For injection powers higher than Hi, the SRL 
exhibits intensity oscillations at approximately the re- 
laxation oscillation frequency wr. For injection powers 
just below SNi, the SRL also exhibits a periodic solu- 



tion which lengthens its oscillation period when raising 
the injection power. This period becomes infinite when 
crossing SNi, where it locks to the injected signal. 

So far standard optical injection behavior has been de- 
scribed. However, the bifurcation diagram in Figure [4] 
reveals the presence of bifurcation lines which are not 
present for a single-mode semiconductor laser [2j. The 
H2 line is a Hopf bifurcation which has a different origin 
than the one encountered in VCSELs [3] and two-color 
lasers [TUj; while the SN2, SN3 and SN4 lines in Figure 
[5] are saddle-node bifurcations, which to the best of our 
knowledge have not been encountered in other optically 
injected laser systems. The H2 Hopf bifurcation is su- 
percritical for values of the detuning close to zero and 
subcritical for more negative values of the detuning (see 
Figure [5| . The stable periodic solution Th associated to 
the supercritical part of H2 disappears when crossing the 
H2-line upwards, where it turns into an injection- locked 
steady state su- This new injection- locked solution Sbi is 
different from s cw which is generated in the saddle-node 
bifurcation. Both solutions are phase-locked to the mas- 
ter laser, but their power distribution amongst the coun- 
terpropagating modes differs; s cw has the optical power 
concentrated in the mode in which we optically inject 
(CW), while su has approximately equal powers in both 
modes. In Figure [5] it can be seen that there is a param- 
eter region in which s cw and su coexist. However, when 
raising the injection power, su disappears at the SN2 line 
while s cw only disappears at the Hi line. For certain de- 
tunings, the H 2 line is located at lower injection power 
than the SNi line, yielding a slightly earlier injection- 
locking (sbi appears at lower injection power than s cw ). 

In the limit cycle Th the phase of the laser field is 
bounded. This means that the phase variables (f>i and 
cj>2 are trapped inside a 27r-wide interval, never crossing 
its boundaries (as opposed to an unbounded or running 
phase solution which freely runs around the phase circle). 
Because the reference frame of our equations is chosen 
to corotate with the phase of the master laser this im- 
plies that although the intensities oscillate, the emitted 
optical frequency is centralized around the master laser 
frequency. The SRL fields are frequency-locked to the 
master laser in Th but not phase-locked. The parameter 
region in which Th exists is partly bounded by H2 where 
it turns into a phase-locked (injection-locked) solution 
Sbi- When moving toward more positive values of the 
detuning Th disappears in a global infinite-period bifur- 
cation (which originates from the TB point, see further 
on) , which is not depicted in Figure [5] When lowering 
the injection power Th disappears in a chaotic attractor, 
approximately at the SN2 line. This transition can be 
seen in Figure [2j There is a clear overlap between the 
region where Th exists and the region confined between 
SNi and Hi where s cw exists, yielding the coexistence of 
an injection-locked solution and a frequency-locked limit 
cycle. 

The SN 2 , SN 3 and SN4 lines are three resonance 
tongues. In the case of an optically injected single mode 
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FIG. 5. (Color online) Bifurcation diagram in the (A, Ei) 
plane. This is a blow-up of Figure[4]in order to show the three 
new saddle-node bifurcations (SN2, SN3 and SN4), which cor- 
respond to three additional resonance tongues. The curves 
SN2 and SN4 have their origin symmetrically around A = 
(originating from the bidirectional solutions), while the SN3 
curve has the same origin as the SNi curve of Figure [4] (origi- 
nating from the unidirectional solutions). The notations s cw , 
s C cw, S),i and Th indicate where these respective solutions 
(co)exist. The shaded area indicates regions of chaotic be- 
havior. Conventions as in Figure [4] 

semiconductor laser we are only faced with a single res- 
onance tongue, corresponding to our SNi line. The oc- 
currence of four different resonance tongues in the case 
of a SRL can be explained by the presence of four diffe- 
rent steady state solutions for the solitary semiconductor 
ring laser at the bias current value we have chosen [T7] , 
In the unidirectional regime the solitary SRL has two 
stable unidirectional solutions and two unstable bidirec- 
tional solutions. The frequencies of these unperturbed 
solutions, which we calculated in Section [IT] can now be 
matched to each resonance tongue since the origin of a 
resonance tongue is located at the frequency of the orig- 
inal non-injected solution. 

For weak injection, the SNi and SN 3 lines are located 
at a detuning A = Acw = Accw (this can be verified 
in Figure [5] see eq. Q for the numerical value), indi- 
cating that these bifurcation curves are related to the 
two unidirectional modes. The SNi line corresponds to 
the CW solution in which we inject and yields the stable 
locking boundary, while the SN3 line corresponds to the 
CCW solution. In the same way a detuning A = Aip 
and A = A p (see eq. ([5|) indicates that the SN 2 and 
SN4 lines respectively correspond to the in-phase and the 
out-of-phase bidirectional solution. Note that these bidi- 
rectional solutions are both unstable for the solitary SRL, 
so for low injection powers the SN 2 and SN4 lines corre- 
spond to bifurcations of unstable structures. Neverthe- 



less part of the SN 2 line is a bifurcation of a stable struc- 
ture, more precisely a bifurcation in which the injection- 
locked solution su disappears, as mentioned before. 

The presence of the SN3 line confirms the intuitive 
reasoning that due to the symmetry of the solitary SRL 
(stable CW and CCW states) and the phase-coupling 
between the fields, two separate injection- locked states 
should originate from the CW and the CCW solution at 
low injection power. When crossing the SN3 line from 
below, an injection-locked solution s ccw appears near the 
original CCW solution through a saddle-node bifurca- 
tion. It is the CCW equivalent of the s cw solution asso- 
ciated to the SNi line, but it has a much smaller basin 
of attraction. The SN 3 line is both steeper and trun- 
cated at the top compared to the SNi line. The increased 
steepness can be understood from the backscattcring phe- 
nomenon. The optically injected light that effectively 
couples into the CCW mode does so through backscat- 
tering. Therefore, the amount of optically injected light 
coupled into the CCW mode is smaller than that coupled 
into the CW mode, increasing the amount of optically 
injected light needed to phase-lock the CCW mode. For 
higher injection powers (higher Ei), s ccw becomes un- 
stable because of the increased effective gain of the CW 
mode when optically injecting it. 

Finally, the SN4 line is a bifurcation of the out-of- 
phase bidirectional solution, which is an inherent unsta- 
ble structure. Crossing the SN4 line does not lead to a 
readily observable change in the system dynamics. 

In Figure [5] we can also see that the Hopf curve H2 and 
the saddle-node curve SN 2 become tangent to each other 
at two different points, so-called codimension-2 bifurca- 
tion points. An analysis of the eigenvalues near these 
points reveals that the TB point above the stable locking 
region is a Takens-Bogdanov point, yielding a double- 
zero bifurcation at that point [59]. It follows from the 
local theory that there must be a global (homoclinic) bi- 
furcation line originating from the TB point somewhere 
between H 2 and the subcritical part of SN 2 . This bifur- 
cation line is not visible in Figure [5] but it is responsible 
for the disappearing of the limit cycle T h , as mentioned 
before. 

The GG point is a Gavrilov-Guckcnhcimcr point, yield- 
ing a fold-Hopf (or zero-pair) bifurcation [29 , only 
involving the SN 2 and the H 2 bifurcation lines (the 
SNi line in Figure [5] only seemingly crosses the GG 
point). The nature of the GG point, which is located at 
Agg ~ —1.645 ns _1 , implies that the H 2 Hopf bifurca- 
tion changes stability when passing GG while increasing 
A [23]. In our case this means that for A < Agg, the 
limit cycle created by (and below) H 2 is unstable, while 
for A > Aqg it is stable. Moreover, the unstable cy- 
cle originates from the unstable node generated by SN 2 
(not from sm, from which the stable limit cycle is gener- 
ated) [5S] . This implies that su does not disappear when 
crossing the H 2 line for A < Agg, but only when cross- 
ing the SN 2 line, which is confirmed by the numerical 
simulations. 
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V. ANTI-PHASE CHAOS 

The occurrence of chaotic regimes in optically injected 
laser systems and photonic integrated circuits is a well- 
known phenomenon, and has been found both experi- 
mentally and by numerical simulation [TJ [51 [HI I5"0"M33| . 
For the optically injected SRL we want to focus on a 
particular chaotic regime (see Figure [3|. It is located 
below the stable locking boundary, and at detunings sig- 
nificantly lower than the relaxation oscillation frequency 
lub,. The parameter region in which this behavior can be 
found is indicated in Figure [5] by the shaded area. This 
parameter region has been constructed by interpolating 
points obtained by numerical simulation of equations (JT|) . 
Both the CW and the CCW mode exhibit purely chaotic 
behavior, but nevertheless act in anti-phase. This is illus- 
trated in Figure [3] by the approximately constant value 
of the total optical power inside the SRL (black line). 

We can distinguish two qualitatively different regimes, 
denoted by X and Y in Figure [3j Regime X is an oscilla- 
tion of the counterpropagating modes at a fundamental 
frequency ujao , making it similar to the alternate oscilla- 
tion regime observed in solitary SRLs (see eq. ^) [2"T] . 
Regime Y is characterized by the suppression of the non- 
injected mode (CCW) by the injected mode (CW) while 
oscillating at the double frequency 2u>ao- The time that 
the SRL resides in cither of these regimes also seems to 
be chaotically distributed. 

Unlike other chaotic regimes in optically injected semi- 
conductor lasers, the anti-phase chaotic regime in an op- 
tically injected SRL shows no involvement of any carrier 
dynamics. The variation in carrier inversion is less than 
0.4 %, which is consistent with the constant total power 
in Figure [3] Hence, the chaos is purely due to the bistable 
character of the SRL and relaxation oscillations do not 
play a role in the onset of chaos as it is the case in other 
optically injected laser systems [SJ 1501455] . 

In order to get more insight in the appearance of 
the strange attractor responsible for the chaotic behav- 
ior, we investigate the orbit in a different phase space. 
Ref- [M] introduces a reduced two-dimensional model 
for the solitary SRL which is valid on timescales slower 
than the relaxation oscillations. It consists of a variable 
€ [— 7r/2, 7t/2], a measure for the distribution of the op- 
tical power among the counterpropagating modes, and a 
phase variable ip € [0,27!"], the phase difference between 
the counterpropagating CW and CCW waves. A value of 
8 = — 7r/2 corresponds to unidirectional CW operation, 
while a value of 9 = tt/2 corresponds to unidirectional 
CCW operation. The validity of this two-dimensional 
phase space has been confirmed by experiments on a soli- 
tary SRL HSJUS]. 

Projecting the five-dimensional chaotic orbit onto this 
two-dimensional (9,ip) phase plane yields a phase por- 
trait which is hard to interpret. However, when viewed 
as a Poincare section, it displays considerable structure. 
Since the phase 4>i of the electric field of the CW mode 
is unbounded in this regime, we construct the Poincare 
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FIG. 6. Poincare section corresponding to an orbit of ap- 
proximately 100 fis (~ 16000 points) in the reduced two- 
dimensional (0, if)) phase space of the solitary SRL, when op- 
erating in the regime depicted in Figure [3] 



section by strobing the system whenever <p± is a multiple 
of 2tt (see Figure [6]). The points now fall on a fractal set, 
which can be interpreted as a cross section of the system's 
strange attractor in the anti-phase chaotic regime. 

Note the topological resemblance of Figure [6] to the 
Poincare section of the chaotic behavior of a periodically 
forced Duffing oscillator (see for example Refs. [331155] ). 
A Duffing oscillator is a nonlinear oscillator with a cubic 
stiffness term, which corresponds to weakly damped mo- 
tion on a double well potential |36J. Hence, its symmetry 
is very similar to the optically injected SRL. The peri- 
odic forcing is a sinusoidal term in the Duffing equation 
which makes the system nonautonomous, corresponding 
to the optical injection in the SRL. Furthermore, the dou- 
ble well potential (with two symmetric stable states) has 
its counterpart in the two stable unidirectional solutions 
of the solitary SRL. Incited by this resemblance in sym- 
metry properties and strange attractor topology, we con- 
jecture that the way the optically injected SRL evolves 
to the anti-phase chaos is topologically identical to the 
way the periodically forced Duffing oscillator evolves to 
chaos. The fact that we observe a period doubling of the 
periodic solution just before the onset of chaos seems to 
endorse this hypothesis [36 . 

Since the periodically forced Duffing oscillator is an 
extensively studied archetypical dynamical system, we 
expect that the analogy with the optically injected SRL 
will allow to interpret and predict dynamical regimes of 
operation which are a consequence of the device symme- 
try. 
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VI. DISCUSSION AND CONCLUDING 
REMARKS 

In this paper, we have theoretically investigated opti- 
cal injection in SRLs. Starting from a single-longitudinal 
mode rate equation model for SRLs, we have used nu- 
merical simulations and a bifurcation analysis to reveal 
all the relevant dynamical regimes that unfold for diffe- 
rent parameter values. We have focused on optical in- 
jection where power is injected in only one of the two 
counterpropagating modes, as this is the case in multiple 
experimental setups. The breaking of the intrinsic Z2- 
symmetry of the SRL due to this optical injection leads 
to remarkable differences in its dynamical behavior com- 
pared to other optically injected laser systems. 

Our bifurcation analysis in Section |IV| showed that 
the behavior of one injection-locked solution is similar 
to the injection-locking of other laser systems. It ex- 
ists in the stable locking region which is bounded by a 
saddle-node (for low injected power) and a Hopf bifurca- 
tion line (for high injected power). This injection-locked 
solution has the modal power concentrated in the CW 
mode in which we inject. Hence, it has a high modal 
power suppression ratio. However, the intrinsic bistabil- 
ity of the SRL leads to three separate parameter regions 
in which this injection- locked state respectively coexists 
with a bidirectional injection-locked state (which has a 
lower suppression ratio), a CCW injection-locked state 



(in which the power is concentrated in the non-injected 
CCW mode) and a frequency-locked limit cycle. This 
frequency-locked limit cycle bifurcates into the bidirec- 
tional injection-locked solution through a Hopf bifurca- 
tion when raising the injected power. 

Our numerical analysis revealed a novel anti-phase 
chaotic regime at a value of the detuning which is much 
lower than the relaxation oscillation frequency. It is dif- 
ferent from chaotic regimes observed in other optically 
injected semiconductor lasers because it does not involve 
any carrier dynamics and the total power emitted by the 
SRL remains constant. Furthermore, parallels to the on- 
set of chaos in the periodically forced Duffing oscillator 
will be the subject of future research. 

ACKNOWLEDGMENTS 

This work has been partially funded by the European 
Union under project IST-2005-34743 (IOLOS) and by the 
Research Foundation-Flanders (FWO). This work was 
supported by the Belgian Science Policy Office under 
grant IAP-VI10. G.V. and S.B. are Postdoctoral Fellows 
and L.G. and W.C. are PhD Fellows of the FWO. Fur- 
thermore, the authors thank Prof. T. Erneux (Universite 
Libre de Bruxelles) and Dr. I. A. Khovanov (University 
of Warwick) for stimulating discussions. The authors also 
wish to thank the anonymous reviewer for his or her con- 
structive comments on the manuscript. 



[1] S. Wieczorek, B. Krauskopf, T. Simpson, and D. Lenstra, 

Phys. Rep. 416, 1 (2005). 
[2] S. Wieczorek, B. Krauskopf, and D. Lenstra, Opt. Com- 

mun. 172 (1999). 
[3] G. H. M. van Tartwijk and D. Lenstra, Quantum Semicl. 

Opt. 7 (1995). 

[4] V. Kovanis, T. Erneux, and A. Gavrielides, Opt. Com- 

mun. 159 (1999). 
[5] Z. G. Pan, S. J. Jiang, M. Dagenais, R. A. Morgan, 

K. Kojima, M. T. Asom, R. E. Leibenguth, G. D. Guth, 

and M. W. Focht, Appl. Phys. Lett. 63, 2999 (1993). 
[6] T. B. Simpson, J. M. Liu, K. F. Huang, and K. Tai, 

Quantum Semicl. Opt. 9, 765 (1997). 
[7] M. K. S. Yeung and S. H. Strogatz, Phys. Rev. E 58 

(1998). 

[8] V. Annovazzi-Lodi, S. Donati, and M. Manna, IEEE J. 

Quantum Elect. 30 (1994). 
[9] I. Gatare, M. Sciamanna, M. Nizette, H. Thienpont, and 

K. Panajotov, Phys. Rev. E 80 (2009). 
[10] S. Osborne, A. Amann, K. Buckley, G. Ryan, S. P. 

Hegarty, G. Huyet, and S. O'Brien, Phys. Rev. A 79 

(2009). 

[11] T. Krauss, P. J. R. Laybourn, and J. Roberts, Electron. 
Lett. 26 (1990). 

[12] M. T. Hill, H. J. S. Dorren, T. de Vries, X. J. M. Leijtens, 
J. H. den Besten, B. Smalbrugge, Y. S. Oei, H. Binsma, 
G. D. Khoe, and M. K. Smit, Nature 432 (2004). 

[13] V. R. Almeida, C. A. Barrios, R. R. Panepucci, M. Lip- 



son, M. A. Foster, D. G. Ouzounov, and A. L. Gaeta, 

Opt. Lett. 29 (2004). 
[14] T. Perez, A. Scire, G. V. der Sande, P. Colet, and C. Mi- 

rasso, Opt. Express 15, 12941 (2007). 
[15] L. Gelens, S. Beri, G. Van der Sande, J. Danckaert, 

N. Calabretta, H. J. S. Dorren, R. Notzel, E. A. J. M. 

Bente, and M. K. Smit, Opt. Express 16 (2008). 
[16] G. H. Yuan, Z. R. Wang, B. Li, M. I. Memon, and S. Y. 

Yu, IEEE J. Sel. Top. Quant. 14 (2008). 
[17] L. Gelens, G. Van der Sande, S. Beri, and J. Danckaert, 

Phys. Rev. E 79 (2009). 
[18] L. Gelens, S. Beri, G. Van der Sande, G. Mezosi, M. Sorel, 

J. Danckaert, and G. Verschaffelt, Phys. Rev. Lett. 102 

(2009). 

[19] S. Beri, L. Gelens, M. Mestre, G. Van der Sande, G. Ver- 
schaffelt, A. Scire, G. Mezosi, M. Sorel, and J. Danckaert, 
Phys. Rev. Lett. 101, 093903 (2008). 

[20] E. J. Doedel, A. R. Champneys, T. F. Fairgrieve, Y. A. 
Kuznetsov, B. Sandstede, and X. Wang, Tech. Rep., Con- 
cordia University (1997). 

[21] M. Sorel, P. J. R. Laybourn, A. Scire, S. Balle, G. Giu- 
liani, R. Miglierina, and S. Donati, Opt. Lett. 27 (2002). 

[22] M. Sorel, G. Giuliani, A. Scire, R. Miglierina, S. Do- 
nati, and P. J. R. Laybourn, IEEE J. Quantum Elect. 39 
(2003). 

[23] C. Born, M. Sorel, and S. Y. Yu, IEEE J. Quantum Elect. 
41 (2005). 

[24] J. Javaloyes and S. Balle, IEEE J. Quantum Elect. 45, 



8 



431 (2009). 

[25] M. San Miguel, Q. Feng, and J. V. Moloney, Phys. Rev. 

A 52, 1728 (1995). 
[26] J. Martin-Regalado, F. Prati, M. San Miguel, and N. B. 

Abraham, IEEE J. Quantum Elect. 33, 765 (1997). 
[27] V. Annovazzi-Lodi, A. Scire, M. Sorel, and S. Donati, 

IEEE J. Quantum Elect. 34 (1998). 
[28] S. Wieczorek, B. Krauskopf, and D. Lenstra, Phys. Rev. 

Lett. 88 (2002). 
[29] Y. A. Kuznetsov, Elements of applied bifurcation theory 

(Springer, 2004), 3rd ed. 
[30] J. Sacher, D. Baums, P. Panknin, W. Elsasser, and E. O. 

Gobel, Phys. Rev. A 45, 1893 (1992). 
[31] T. B. Simpson, J. M. Liu, A. Gavrielides, V. Kovanis, 



and P. M. Alsing, Phys. Rev. A 51, 4181 (1995). 

[32] V. Kovanis, A. Gavrielides, T. B. Simpson, and J. M. 
Liu, Appl. Phys. Lett. 67, 2780 (1995). 

[33] M. Yousefi, Y. Barbarin, S. Beri, E. A. J. M. Bente, M. K. 
Smit, and R. Notzel, Phys. Rev. Lett. 98 (2007). 

[34] G. Van der Sande, L. Gelens, P. Tassin, A. Scire, and 
J. Danckaert, J. Phys. B-At. Mol. Opt. 41 (2008). 

[35] S. H. Strogatz, Nonlinear Dynamics and Chaos: with Ap- 
plications in Physics, Biology, Chemistry and Engineer- 
ing (Addison- Wesley, Reading, MA, 1994). 

[36] J. Guckenheimer and P. Holmes, Nonlinear oscilla- 
tions, dynamical systems, and bifurcations of vector fields 
(Springer, 1990), 3rd ed. 



